The role of /-electrons in the Fermi surface of the heavy fermion superconductor /?-YbAlB4 
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We present a detailed quantum oscillation study of the Fermi surface of the recently discovered Yb-based 
heavy fermion superconductor /3-YbAlB4. We compare the data, obtained at fields from 10 to 45 Tesla, to band 
structure calculations performed using the local density approximation. Analysis of the data suggests that /- 
holes participate in the Fermi surface up to the highest magnetic fields studied. We comment on the significance 
of these findings for the unconventional superconducting properties of this material. 

PACS numbers: 71.18.+y, 71.20.Eh, 71.27.+a, 75.30.Mb, 75.10.Lp 



The recent demonstration of superconductivity in the heavy 
fermion compound /3-YbAlB4 Q provided a breakthrough in 
the search for a ytterbium-based analogue of the numerous 
cerium heavy fermion superconductors lllElBll- Theoretical 
models of / electron metals often describe the low tempera- 
ture dynamics in terms of the dense Kondo lattice hamiltonian 
ISll . In this framework, it is expected that hole and electron 
like materials should realize similar ground states. While Yb 
in a 4/^^ electronic configuration is the hole counterpart to 
4/^ Ce, superconductivity in Yb-based materials |6, 7| has 
proved strangely elusive in contrast to those containing Ce, 
until the discovery of superconductivity in /3-YbAlB4. 

Single crystals of ;^-YbAlB4 grow in a plate-like morphol- 
ogy, with a crystal structure corresponding to that of ThMoB4 
ISl . In contrast to the more common square Yb superlattice, 
the Yb ions in /3-YbAlB4 form a lower symmetry distorted 
hexagon, with Yb and Al atoms in the basal ab-plane sand- 
wiched between nets containing an unusual arrangement of 
heptagonal and pentagonal B rings. Alloying with smaller in- 
terstitial atoms favours a smaller 4/^^ ionic configuration with 
the shortest Yb-Yb bond length lying along the c-axis. Ex- 
perimental studies of the magnetic susceptibility |1| show a 
large Ising-like moment along the easy [001] axis, these mea- 
surements demonstrate the presence of local moments at high 
temperatures. 

The occurrence of superconductivity with Tc = 80 mK in 
the cleanest samples of /3-YbAlB4 with po ^ 0.8//f]cm raises 
the intriguing question of how superconductivity develops out 
of this local moment system. Further evidence of novel corre- 
lated electron behaviour (H is signaled by low temperature 
non-Fermi liquid (NFL) behaviour, with the magnetic sus- 
ceptibility diverging for applied magnetic field B \\ c, ac- 
companied by an electrical resistivity p ex T^/^. Applica- 
tion of a small magnetic field at low temperature recovers FL 
behaviour with enhanced specific heat coefficient 7 (gener- 
ally referred to as a heavy FL), suggesting that the /-electron 
bands hybridize with the conduction electrons and thereby 



contribute to the Fermi volume |9|. Individual local moment 
behaviour is found to be recovered above an energy scale 
of T* ~ 200 K 1 1 1. This thermodynamic and transport evi- 
dence for strongly correlated electron behaviour in /3-YbAlB4 
provides impetus to study the electronic structure. 

In this Letter we present a study of the Fermi surface (FS) 
of /3-YbAlB4, observed by measuring quantum oscillations 
in a range of magnetic fields. High quality superconducting 
single crystals with typical dimension 2 x 0.1 x 0.01 mm^ 
were grown according to the method described by Macaluso 
et.al. |8|. Where necessary, electrical contacts were made to 
the crystals using Dupont 6838 silver epoxy, yielding contact 
resistances of ~ 0.25 1^ at room temperature. 

A study of the temperature, field and angular dependence of 
quantum oscillations was carried out using a number of tech- 
niques. Low field Shubnikov-de Haas (SdH) measurements 
were carried out in Cambridge between 10 and 18 T at tem- 
peratures as low as 10 mK, using excitation currents 40 /iA 
below 100 mK to avoid sample heating. The angular depen- 
dence from B II [001] to [100] and B \\ [001] to [010] was 
studied at low fields. High field measurements were carried 
out using the 45 T hybrid magnet facility at the National High 
Magnetic Field Laboratory in Tallahassee, where oscillatory 
signals were observed in the field range 35 - 45 T, to temper- 
atures as low as 40 mK. Rotations from 5 || [001] to [100] 
were studied using a capacitance based torque magnetometer 
to observe the de Haas-Van Alphen (dHvA) effect, while ro- 
tations from B II [001] to [010] were observed using the SdH 
effect. The DC magnetization was measured using a pulsed 
field at the ISSP in Tokyo. 

Fig. [T] shows sample oscillations in raw data treated with 
a polynomial background subtraction, the inset of the upper 
(lower) panel corresponds to high (low) field measurements 
with B II [010]. The main panels show Fourier transforms 
of the data at various temperatures, with several peaks clearly 
visible in the spectrum associated with extremal orbits of the 
FS. Estimates of the quasiparticle mean free path I free were 
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FIG. 1 : (color on[inc)Inset: Percentage change in the sample resis- 
tance as a function of fiQld.Top panel: SdH data obtained from 35 - 
45 T B|| [010] and associated 1/B FFT spectra at various tempera- 
tures. Lower panel: SdH data obtained from 15 - 18 T B|| [010] and 
associated 1/B FFT spectra at various temperatures. 

obtained from field dependent amplitude fits to the Dingle fac- 
tor for impurity scattering Rd = exp(Cr/2//^ee) (where Cr 
is the circumference of the associated real space orbit). Fit- 
ting yield I free = 3400A for orbits in the basal plane and 
I free = 460 A for orbits perpendicular to [100], assuming 
circular orbits. The mean free path is of comparable length 
to that calculated by more indirect means |10|, placing (3- 
YbAlB4 in the clean limit for superconductivity. 

In order to clearly elucidate the role of the /-electrons in (3- 
YbAlB4 experimentally measured quantum oscillations were 
compared with predictions from band theory for the 4/^^ case 
of /3-LuAlB4, where the /-shell is filled and falls below the 
Fermi level, and for the 4/^^ case of /3-YbAlB4, where the /- 
hole is believed to be itinerant. Ab initio electronic structure 
calculations were performed using the Wien2K package ifTTIl . 
Calculations were conducted under the local density approxi- 
mation including the effect of spin-orbit coupling (LDA+SO) 
using full potentials and a LAPW basis. A grid of 39 x 39 x 64 
k-points was used for the Brillouin zone integration. For com- 
pleteness a ferromagnetic LDA+SO+U calculation with the 
moment along c-axis has also been performed for /3-YbAlB4, 
yielding essentially the same FS as that found for /3-LuAlB4. 

Fig. [2] shows a comparison of the angular dependence of 
the oscillation frequency of our experimental data at high 
and low fields plotted alongside the angular dependence of 
the LDA+SO FS for p-LuA\B^ (top panel) and /3-YbAlB4 
(bottom panel). Fig. [3] shows the calculated FS for the main 
bands that cross the Fermi level for both the /?-LuAlB4 and P- 
YbAlB4 cases. A strongly 3D FS is predicted by the LDA+SO 
calculations despite the apparent 2D nature of the crystal 
structure, due to strong orbital hybridization in the short c- 
axis direction of the unit cell. 

The lower panel of Fig. |2] corresponds to the itinerant hole 
case, and shows a more satisfactory match between experi- 
ment and theory. In particular there is reasonable agreement 



1 1 1 1 


1 1 1 1 

LUAIB4 

□ □ I 




^ — ""^ ^ ^ 

— : 

□ 

^^^^ ES^^ 



1 1 1 1 


1 1 1 1 

YbAIB4 




□ □ I 

iwfr^^ ^^^^ □ 


- a ' — 





-80 -60 -40 -20 20 40 60 80 
[100] [001] [010] 



Angle (Degrees) 

FIG. 2: (color online) Angular dependence of experimental quantum 
oscillation frequencies at low (filled black square) and high (unfilled 
black squares) field compared to that of the LDA+SO FS: 
Top panel: /3-LuAlB4. Green^band 91 and purple^band 93 
Bottom panel: /3-YbAlB4. Red^band 89 and blue^band 91 



TABLE I: Experimentally observed orbits and masses determined by 
fitting to the Lifshitz-Kosevich relation for a single mass. Together 
with the band mass of the assigned orbit where a corresponds to 
band 89 and /3 to band 91 where an assignment was possible, 7 cor- 
responds to experimental frequencies that are not assigned. 
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with the experimental frequencies observed with B \\ [001], as 
well as for the frequencies below 500 T observed throughout 
the rotation. The inclusion of the spin-orbit interaction was es- 
sential for obtaining quantitative agreement along the c-axis, 
which is indicative of the importance of the 4/ moment at the 
Fermi level. The large frequencies predicted for 5 || [010] and 
B II [100] are greater than those observed experimentally, yet 
possess similar angular dependence. The discrepancies are 
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FIG. 3: (color online) Panels a) and b) show the FS corresponding 
to bands 91 and 93 in /3-LuAlB4 (with a filled / shell), calculated 
using the LDA technique described in the text. Panels c) and d) show 
analogous calculations for bands 89 and 91 in ;5-YbAlB4, assuming 
itinerant /-electrons. In all cases effects of spin-orbit coupling were 
included. Quantum oscillation frequencies correspond to the area 
enclosed by extremal orbits, for example we ascribe f3i to the belly 
of band 91 around the Z point. Two additional small pockets from 
band 93 of ;^-YbAlB4in the plane of Z are not shown in the figure. 

likely indicative of the presence of correlation effects beyond 
the LDA. The high field data, taken up to 45 T, shows broad 
agreement with the angular dependence of the low field data. 

The full bandstructure resulting from the LDA-fSO calcula- 
tion of /3-YbAlB4 along with the associated density of states 
is depicted in Fig. |4] Spin-orbit coupling splits the / level 
into J = 7/2 and J = 5/2 components separated by 1.4 eV, 
with crystal electric field splitting further lowering the sym- 
metry of the ground ^^7/2 state, thereby yielding a ground 
state doublet. By summing over the projected density of states 
for J = 7/2, an occupation of 7.42 was obtained, implying an 
/-hole value of 0.58. This mixed valence character yielded by 
the LDA+SO solution is inconsistent with the experimentally 
observed local moment, implying that strong coupling effects 
beyond the LDA are important. 

Visualization of the charge density from the itinerant calcu- 
lation indicates that the predominant bonding charge density 
lies about the Boron layer. The projection of the Yb J = 7/2 
density of states from the full density in Fig. |4] demonstrates 
the strong hybridization of the Yb / level with states from 
other atoms, in this case the Boron layer, this is of importance 
to the itinerant electron character of this material. 
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FIG. 4: Main panel (color online): Band structure produced by 
LDA-i-SO. The inset shows the Wigner-Seitz cell used to define the 
crystallographic directions. Right panel: The associated density of 
states (DOS), with the total DOS in black and the contribution from 
the J = 7/2 states in red. 



We proceed to compare the measured quasiparticle effec- 
tive mass (m*), extracted by fitting the temperature depen- 
dence of the FFT spectral weight to the standard Lifshitz- 
Kosevich relation (Fig. [5]), with the calculated band mass. Ta- 
ble 1 shows a comparison of experimentally measured masses 
at low fields and their enhancement over calculated band 
masses. The sizable enhancement in effective mass is con- 
sistent with significant /-electron hybridization at the Fermi 
energy. An estimate of the contribution to the electronic spe- 
cific heat component from the experimentally observed FS 
sheets is obtained by applying the experimental mass renor- 
malization to the density of states calculated for the observed 
sheets. While the total density of states at the Fermi energy 
from electronic structure calculations is small, equivalent to 
an electronic specific heat component of 7 = 6.4 mJ/mol K^, 
taking the average of the experimental enhancement yields 
7 ~ 50 mJ/mol - comparable to the experimentally de- 
termined value of 7 ~ 75 mJ/mol measured at 9 T. As 
such we are able to account for the enhancement in the spe- 
cific heat, furthermore this suggests that there is no large dif- 
ference in effective mass between spin channels. 

The large density of states associated with hybridized bands 
in /3-YbAlB4 raises the possibility that the paramagnetic HFL 
can be polarized at modest applied magnetic fields. Polar- 
ization of the HFL would in turn lead to a shift in chemical 
potential, associated with a reduction in the experimentally 
measured effective mass. We therefore further probe the ex- 
tent of /-electron involvement in the FS up to high magnetic 
fields. Fig. [5] shows the effective mass evolution as a func- 
tion of the applied magnetic field. Rather surprisingly, the 
modest mass enhancement persists at high magnetic fields, in 
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FIG. 5: Top panel: Magnetization measured at 1.4 K B || [001] up 
to 68 T together with the field dependence of the 04 orbit B \\ [010]. 
Note a superlinear increase in M over the range quantum oscillations 
were observed. Bottom panel: Lifshitz-Kosevich fit to fS^ orbit at 
15.7 T for a single mass. 



fact the effective mass is seen to increase with magnetic field 
measured for 5 || [010]. From the energy scale E = /ieff^ 
associated with the highest applied magnetic fields at which 
quantum oscillations have been observed, we infer that /- 
electrons are involved in the Fermi surface at least up to an 
energy scales 100 K. 

In order to understand the increase in effective mass with 
magnetic field, we consider the trend in magnetization at these 
fields (Fig.[5|. Measurements of the magnetization up to 68 T 
reveal an absence of a metamagnetic transition 1121 [TSl . The 
subtle non-linear increase in magnetization at higher fields 
implies a non-linear change in the density of states around 
the Fermi level. Accordingly, the enhancement of the effec- 
tive mass shown in Fig.|5]is consistent with an increase in the 
magnetic susceptibility. 

Our mapping out of the electronic structure in /3-YbAlB4 
could provide useful insight to the origin of unconventional 
superconductivity in this material. Comparing the /?-YbAlB4 
system with the prototypical Ce-based unconventional super- 
conductor CeCoIns [3, JL4|, we note potential similarities. 
Both these systems are located in proximity to putative quan- 
tum critical points in the vicinity of which NFL behaviour is 
seen; while CeCoIns lies on the antiferromagnetic side, /3- 
YbAlB4 is positioned on the paramagnetic side. Our exper- 
imental measurements in /3-YbAlB4 are performed at high 
magnetic fields beyond the region of NFL behaviour, enabling 
us to access the low energy excitations of the FL prior to its 
evolution to the quantum critical regime. We find direct ev- 
idence that /-electrons participate in the Fermi surface of P- 
YbAlB4 right up to temperatures of the order of 30 K, reflect- 
ing the large Kondo lattice energy scale in this material Q - 
in fact similar to CeCoIns 115 J . 



Given the likely association of an enhanced coherence en- 
ergy scale with the appearance of superconductivity in these 
materials, it appears unexpected that the superconducting tem- 
perature in /3-YbAlB4 is considerably lower when compared 
to the Kondo temperature scale than in the case of CeCoIns. 
Potential clues as to the suppressed Tsc in /3-YbAlB4 arise 
from the FS topology. An outstanding question relating to 
the appearance of unconventional superconductivity in heavy 
fermion systems relates to the origin of Cooper pairing from 
within these large moment Kondo lattice systems. Despite 
the large local moment in /3-YbAlB4, our quantum oscillation 
measurements find FS pockets in which itinerant /-electrons 
are involved. We are thereby able to pinpoint the locations 
in reciprocal space comprising heavy electrons from which 
superconducting Cooper pairs likely originate. Among rea- 
sons for the vastly lower superconducting temperature in (3- 
YbAlB4 may perhaps be the robust 3D nature of these FS 
pockets: it has been suggested that the size of Tsc is enhanced 
in 2D |16|. 

An interesting question for further exploration involves the 
nature of evolution of the observed heavy FS pockets as the 
system is tuned toward the other side of the quantum critical 
point, where magnetism may potentially replace superconduc- 
tivity as the groundstate. We find that an antiferromagnetic 
LDA+U calculation in fact yields a lower energy ground state, 
hinting that the system lies close to an antiferromagnetic insta- 
bility. The resolution of these questions and how /?-YbAlB4 
fits into the Kondo lattice phenomenology is anticipated to 
form the basis of future investigations. 
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